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Abstract 

We discuss in this paper a method of finding skyline or non-dominated points in a set P of np 
points with respect to a set S of ng sites. A point Pi G P is non-dominated if and only if for each 
Pj £ P, j ^ i, there exists at least one point s £ S that is closer to Pi than pj. We reduce this problem 
of determining non-dominated points to the problem of finding sites that have non-empty cells in an 
additive Voronoi diagram with a convex distance function. The weights of the additive Voronoi diagram 
are derived from the co-ordinates of the points of P and the convex distance function is derived from 
S. In the 2-dimensional plane, this reduction gives a 0{{ns + np) logng -I- np lognp)-time randomized 
incremental algorithm to find the non-dominated points. 

1 Introduction 

Consider a trip to a conference in a new city! A set P of np hotels (located at fixed locations) has already 
been identified from a travel guide. On reaching the city, the scientist identifies a set S of ns sites to visit, 
say for example the conference venue, museum, restaurant, garden, beach, etc. The scientist wants to visit 
all sites in S but prefers a hotel that has at least one site in S that is closer to it than any other hotel. Now, 
which are the most interesting hotels in the set P with respect to the sites of S in terms of distance? A 
hotel is interesting if it has at least one site closer to it than any other hotel. This problem gives rise to the 
spatial skyline queries [11 . A point G P is a skyline point if it has at least one site in 5* that is closer to 
Pi than to any other point in P. 

There can be applications of this problem in other areas like identifying a set of buildings for quick evacuation 
in case of multiple fires. Here, the set of buildings is P and the set of multiple fires is S. The set of skyline 
points are the buildings among P that are to be evacuated ahead of the other buildings. Sharifzadeh and 
Shahabi |llj identify some other applications as well. 
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1.1 Formal Definition 



Let P — {pi, . . . ,Pnp} be a set of np points and S = {si, . . . , Snsl be a set of ns sites in M^. Let d{x, y) be 
the usual L2 distance in IV^ . Each point pi £ P has ns spatial attributes determined by all the ns distances 
d-iPij e S), of the points in 5" to p,;. Next, we define domination. 

Definition 1 (Domination) 11 If Given a set P of np points and a query set S of ns points in the plane, 
Pi (z P dominates pj G P (j ^ i) with respect to S if and only if d{pi, s) < d{pj, s), Vs G 5*. 

In this setting, the distances to the points of S can be considered feature vectors that describe the sites Pi 
and Pj . The feature vector for pi dominates the vector for pj if and only if it dominates on all coordinates. 
li Pi dominates pj, pj is a non-interesting point (hotel) vis-a-vis pi with respect to S. Note that, pi is not 
dominated by pj if it has at least one point in S that is closer to it than pj . 

We define skyline points as the set of those points in P which are not dominated by any other point in P 
with respect to S. We denote the set of skyline points as SV. We will use the terms skyline points and 
non-dominated points interchangeably. 

Definition 2 (Skyline point) pi € P is a skyline point, if and only if we have the following: 

for each pj (z P, i ^ j, 3s G S such that d{pi, s) < (pj, s) (1) 



In our problem, we want to extract skyline points of P with respect to S. Consider a brute force approach. 
Let h{pi,pj) be the half-plane of pi with respect to pj. For each pi, determine if there is at least a point 
s G 5 which lies in h(pi,pj) for all j =/= i. If at least one s is found for every pj, then pi is a skyline point. 
This takes 0{npns) time for each pi. With the assumption that np = ns ~ n, the total time complexity is 
0(n3). 

1.2 Prior Work 

Given two points pi — {pj,pf,... ,pf) and pj — . . . ,Pj) in Ff^, pi dominates pj if and only if < 

for 1 < m < d and < fo^' some 1 < m < d. For a point set P in M'^ the skyline query finds those 
points in P which are not dominated by any other point. Skyline operator was introduced by Borzsonyi et al. 
[1]. They implemented skyline query to update an existing (relational, object-oriented or object-relational) 
database system with a new logical operator that they refer to as the skyline operator. Borzsonyi et al. 
used divide-and-conquer techniques and index structures to solve the problem in 0(nlog''^^ n + nlogn) time 
where d is the number of dimensions of the points. Since the introduction of skyline query and skyline 
operators by Borzsonyi et al. ^ , there have been several works using nearest neighbor search [7] , sorting ^ 
and index structures [9l [13] . These works mostly try to show an improvement over the results of Borzsonyi 
et al. \X\ experimentally. 

The problem of Spatial skyline query, as introduced in this paper, was first addressed by Sharifzadeh et al. 
[TT] . The basic difference between the work of Borzsonyi et al. pLj and Sharifzadeh et al. [11] lies in the 
definition of domination. Borzsonyi et al. define domination between two points based on their respective 
coordinates, whereas, Sharifzadeh et al. define domination between two points with respect to a set of points 
as given in Section [O] Note that, the method of Borzsonyi et al. can be applied to the problem of spatial 
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skyline query, but then the time complexity would be 0(np log"^~^ rip + rip log np) if all the 0{npns) 
distances are already computed. Sharifzadeh et al. [llj propose an 0(n|.|iS7^| + ^/np) algorithm for the 
above problem where \S'P\ denotes the cardinality of the solution set. They solve this problem using Voronoi 
diagram, convex hull and delaunay graph. Observe that, if we put ns = 0{\) in the time complexity derived 
by Sharifzadeh et al. [11], their worst-case time complexity becomes 0{np). Moreover, Son et al. [12] 
have shown that the algorithm and the time complexity analysis of Sharifzadeh et al. jllj is incorrect and 
proposed a solution whose time complexity is 0{np\SV\\og \C1-L{S)\ + nplognp) where ISV] and \CH{S)\ 
denote the cardinality of non-dominated points and convex hull of the set S respectively. If we look at the 
worst case complexity of the algorithm devised by Son et al. [T2] , it turns out to be 0{np log ns + np log np) . 



1.3 Our Work 

As discussed earlier in Section [l~T| the set of skyline points in P is the non-dominated subset of P with respect 
to S. We first show using lifting techniques [2] that the set of non-dominated points has a correspondence 
with the lower envelope of cones. The cones are obtained as follows. Choose an origin in the convex hull of 
S. Consider an unit paraboloid with its apex at this origin. Fix a point p € P. Let its lifted version on the 
unit paraboloid be p' . Now, for all points s e CH{S), consider the discs centered at s and passing through p. 
Each such disc, if lifted onto the unit paraboloid, forms a plane passing through p'. Now, the lifted versions 
of all such discs form a cone with apex at p' . So, we will have np such cones. In Section [2.1| we show that 
the non-dominated points of P correspond to the apices of the lower envelope of the said cones. 

This scheme can be alternatively interpreted in a Voronoi diagram model as follows. The lifted coordinates of 
the points in P are taken as their additive weights. With the same origin and the unit paraboloid as before, 
lift each point of S to planes tangent to the unit paraboloid. Translate these planes to include the origin and 
intersect their halfspaces to define a cone. The cross section of this cone at a unit distance vertically above 
the origin defines a convex polygon. We take this convex polygon to define a convex distance function. Now, 
the lower envelope of cones correspond to an additively weighted Voronoi diagram for the convex distance 



function. We show in Section 2.2 that the skyline points of P with respect to S are those with non-empty 



Voronoi cells under this convex distance function determined by S with additive weights determined by P. 

After having shown the relation of skyline points to the non-empty cells of the said Voronoi diagram, we 
proceed in Section [s] along the lines of McAllister et al. [H] , where the authors have shown that compact 
diagrams that avoid high combinatorial complexity of Voronoi diagrams under convex distance function can 
be used for solving certain problems. We show that computing such a compact diagram which can be used 
to find non-empty Voronoi cells under a convex distance function determined by S with additive weights 
determined by P takes 0{{ns + np) logng + np log np) time for a randomized incremental construction. 



2 Reduction to a Voronoi diagram 

In this section we relate point domination to additively-weighted Voronoi diagrams of a convex distance 
function in IR^. We will define these terms as we go, culminating in the following: 

Result 1 The non-dominated points of set P with respect to sites S are those with non-empty Voronoi cells 
under a convex distance function determined by S with additive weights determined by P. 

A brief sketch relating non-domination to non-empty Voronoi cells under a convex distance function with 
additive weights is as follows. Choose an origin in the convex hull of S. Lift each point of S to planes 
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tangent to the unit paraboloid and sites P to points on this paraboloid and take the lifted coordinate 
as an additive weight. Translate the planes to include the origin and intersect their halfspaces to define a 
cone, from which we take the convex distance function. Thus, the convex distance function is determined by 
S. The lower envelope of cones with apices at the lifted points of P bounds the non-dominated (additively 
weighted) points, and can be interpreted as an additively weighted Voronoi diagram for the convex distance 
function. We elucidiate further. 

2.1 Dominated points and the cone 

Let C{x,y) denote a disc with center x and radius equal to d{x,y). For each point pi G P, consider discs 
C{sk,Pi) centered at each G 5*. Obviously, for any point p £ inside the disc C{sk,Pi) Sk is closer to p 



than to Pi. Therefore, if a point pj dominates pi, then pj £ || C{sk,Pi)- Let Di — || C{sk,Pi). We term 



Di as the dominator region of pi. The significance of such a dominator region is that any point belonging to 
Di dominates pi with respect to S. Now we have the following observations. 

Observation 1 For any point pi £ P: 

(i) Di is nonempty and may overlap with another Dj where j ^ i 

(ii) Di is a convex region bounded by circular arcs. 

(iii) Pi is a skyline point if its dominator region Di does not contain any point pj £ P,j ^ i in its interior. 

(iv) ij a point Pj £ P lies inside the region Di,i ^ j, then the dominator region Dj of pj is a subset oj Di, 



Proof 1 (i), (ii) and (iii) are trivial. We prove (iv) by contradiction. If Pj lies inside Di then all Sks' 
are nearer to Pj than pi; so pj dominates pi with respect to S. Assume that Di does not contain Dj. This 
implies that there is a point z £ Dj and z ^ Di, such that z dominates pj but not pi. But pj dominates pi. 
As the domination relation is transitive there is a contradiction. Hence, we have the observation. 

Note that, for any dominator region Di of a point pi, the boundary of Di is determined by at most ns 
circular arcs. So the total complexity of this configuration of dominator regions for all the points in P can 
be O^ngnp) ~ O(n^) under the assumption that ng = np ~ n. 

Let C'H{S) denote the convex hull of S. Assume that the origin for the point sets S and P lies inside C'H{S). 
Now, lift each point of S and P to points on a unit paraboloid where {'i' = {x,y,z)\z ~ x'^ + i/^}[2]- 
So, any point p = {x,y) in the plane is lifted to a point p' = {x,y,x'^ + y^) on ^. We will refer to this 
geometric transformation as paraboloid or lifting transformation. For any and Pi in the plane, let s'f, and p[ 
respectively denote the lifted image on the unit paraboloid ^. Now, a paraboloid transformation of a circle 
C — {{x—ci)^ + {y—C2)^ = r^} in the two dimensional 2:y-plane is a curve C" = {z—2cix—2c2y+c\+C2 — r'^ — 
0} on ^ [To]. The equation of C depicts a plane in 3-dimensional space. Observe that the disc C{sk,Pi) is 
a plane in 3-d; we denote this plane as C'{sk,Pi). Moreover, each C'{sk,Pi) (fc = 1, . . . ,ns) passes through 
p'^ and is parallel to the tangent plane of the unit paraboloid ^! at sj.. Now we have an observation linking 
the dominator region Di of pi with a cone having its apex at p^. 

Observation 2 For any point pi, the dominator region Di in the 2D plane is mapped to a cone fi^ in S-space 
with its apex at Pi under paraboloid transformation. 





i.e. Dj C Di 
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Proof 2 From Observation 1 

through p[ where C'{sk,Pi) and p'^ are the lifted images of C[sk,Pi) and Pi respectively by paraboloid trans- 
form,ation. The intersection of the upper half spaces o/ {C"(sfc,pi)| \/k}'s defines a cone Oj with apex at p[ 
in 3-space. This cone in 3-space corresponds to the dominator region in the plane. 

Let pj dominate pi. We already know from Observation [l] that Dj C Di. Next, we explore the relation 
between fli and flj where flj is the cone corresponding to the dominator region Dj of pj . 

Lemma 1 The cone flj corresponding to the dominator region Dj of point pj will contain the cone fli 
corresponding to the dominator region Di of point pi, if Pj dominates pi.. 

Proof 3 Consider a point s'^ on the unit paraboloid corresponding to s„i G S . Tm is the tangent plane 
of at sj„. Translate the plane t,„ to include p'^ and we denote that plane as r™. The projection of the 
intersection of t™' and ^' on the plane is C{srn,Pi)- As pj dominates pi, during translation t„i will include 
Pj before p^. This is true for all Tmijn = 1, . . . ,ns)- Since f2i is defined by intersection of the upper half 
spaces of {r™| Vm = 1,2,..., ns}, fij will contain fli. 

Observe that a point can be dominated by more than one point. From Lemma [T] we can conclude that a 
cone f2i in 3-space corresponding to a dominated point pi will be contained in atleast one cone flj . We could 
relate till now the relation between containment of cones and domination between two points. An extension 
encompassing all points in P indicates the relation between non-domination, and hence, skyline points, to 
lower envelopes of cones. Now we state the following theorem. 

Theorem 2 The apices corresponding to the lower envelope of the cones {fii, . . . ,il„p} are the nondomi- 
nated points, i.e the skyline points of the point set P. 

Proof 4 Follows from Observations^ and^ and Lemma\^ 

2.2 Relation of lower envelope of cones to additively weighted Voronoi diagrams 
of a convex distance function 

As deduced in Theorem [2] the set of skyline points is nothing but the points of P corresponding to the apices 
of the lower envelope of the cones fi^. Constructing this lower envelope of cones is costly. So, we try to explore 
a relation between additively weighted Voronoi diagrams of convex distance functions with lower envelope of 
such cones. To give the details, we need to define convex distance functions and additively-weighted Voronoi 
diagrams. 

Definition 3 (Convex distance function) Minkowski showed that any convex set M whose interior contains 
the origin defines a convex distance function dM{p,<l), where the distance from point p to q with respect to 
M is the amount that M must be scaled to include q — p. 

d-AiiP, q) = inf{A >0:q-peXM} 



know that Di = C{sk,Pi)- Now each C'{sk,Pi) is a plane passing 



fc=i 
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If M is closed, then the infimum operation can be replaced by the minimum operation. A distance function 
may not be a metric, since (Im is symmetric if and only if M is centrally symmetric. (If we denote the 
reflection of the set M through the origin by Af®, then dM{p,q) = dMe{q,p).) The distance function (Im 
does always satisfy the triangle inequality for points [5]: dM{p,q) + dMiq,r) > dM{p,r). Note that the 
boundary of M serves as the unit ball for this distance function. For a fixed p, the graph of dM{p, q) as a 
function of g is a cone with apex at p. 

Definition 4 (Additively-weighted Voronoi diagram) Given a finite set of sites P = pi,p2, • ■ ■ ^Pn C M'^ , 
with additive weights lui,u}2, ■ ■ ■ t^u, o,nd any distance function, d{p,q), we can define a generalized Voronoi 
diagram by labeling each site q G M'^ with its set of closest sites, 



and partitioning the plane into maximally connected regions having the same labels. Voronoi cells are regions 
with a single closest neighbor and vertices have degree d+1 (or more, in degenerate configurations). 



Figure 1: A Voronoi diagram of 6 sites in the plane using a convex quadrilateral as distance function; its 
view as a lower envelope of cones. 

Figure [ija) illustrates a simple example for 6 distinct sites, all having weight zero. The distance function 
is the black convex quadrilateral around the point at the origin, and each cell is drawn in a different color. 
Figure [ijb) shows the cones for which the Voronoi diagram is the lower envelope. When all weights are zero, 
each distinct site has a non-empty cell. The following observation, which we state without proof, is central 
to our idea of relating the lower envelope of cones with additively weighted Voronoi diagrams of a convex 
distance function. 

Observation 3 The cell for a site shrinks if we increase the weight of the site; we essentially translate the 
cone upwards until the lifted site disappears from the lower envelope. 



label{q) — argmind{pi,q) + uj., 

iG[l,n] 




(a) 
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2.3 Reduction 



We can do the reduction by simply giving a different interpretation of the inequahties defining dominance. 
Recall that, as per Definition [l| a point p dominates q {p,q £ P) with respect to the set of sites S if and 
only if Vs g S, d{p, s) < d{q, s). 

Without loss of generality, choose the origin to be some point inside CTL{S), then assign each site p G P a 
weight bjp = p.x^ +p.y'^ , which we can use as the lifting coordinate. If we square both sides of the dominance 
inequality we obtain a linear expression in site coordinates and weights: 

(p - s) • (p - s) < {q-s)-{q~s) (2) 
if and only if ujp < ujq — 2s ■ (q — p) . (3) 

Note that, for each s G S*, the above inequality gives rise to a (hyper)plane passing through p' (= (p.x,p.y, Wp)), 
which is the lifted version of p. The intersection of the halfspaces defined by the points of S gives a cone 
with apex p' that contains all points dominated by p (see Lemma [T]). We now show that the collection of all 
cones for the points P corresponds to an additively weighted Voronoi diagram for a convex distance function 
defined by S. 

We claim that if we intersect this cone with a (hyper)plane u> = Up + 1, then we obtain a convex polytope 
that defines a distance function containing (p.x^p.y,ujp + 1) as its origin. Let M be the projection of this 
convex polytope onto the (hyper)plane a; = 0. Note that, moving from this origin ((p.x,p.y,U!p + 1)) by 
any vector (v.x,v.y^O) must leave the polytope, since any v can be expressed as a convex combination 
V — —'^^i<:i<:ns '^i^i where reals a-i are not all equal to zero because the points S = ni, n2, . . . ,ns contain 
the origin in their convex hull. Thus, the Voronoi diagram of distance dM for sites p E P with weights ujp 
generates the same lower envelope of cones. The sites that are not dominated with respect to S are those 
with non-empty cells. From the above discussion, and Lemma [l] and Observation [3] we get the final result 
that was stated in Result [T] 

Theorem 3 The skyline or non- dominated points of set P with respect to sites S are those with non-empty 
Voronoi cells under a convex distance function determined by S with additive weights determined by P. 

3 Computing non-dominance in the plane 

Theorem[3]tells us points in P that have non-empty Voronoi cells in an additively weighted Voronoi diagram of 
P under a convex distance function determined by S correspond to the skyline points. We have also described 
the method to obtain this convex set M from S at the beginning of Section [I Note that |A/| = 0(|5|). 
So, our problem of finding skyline points is now transformed into computing non-empty Voronoi cells in 
an additively weighted Voronoi diagram under a convex distance function M. Actually, we want to avoid 
computing all the edges because, as seen in Figure [ij a), they are fairly complicated, even when all weights 
are zero. The bisector between two sites can consist of 0(|M|) line segments, so the generalized Voronoi 
diagram of sites P can then have complexity 6(|-P| • \M\). We can see this from the cone view. So, we use the 
compact piece-wise linear Voronoi diagram concept of McAllister et al. [H]. A compact Voronoi diagram is 
an approximated version of an abstract Voronoi diagram(AMT)) ffT defined using a convex distance function. 
An AVD is defined only in terms of bisectors of pairs of sites and are computed using the ordering of the 
two points along a bisector and the ordering of these bisectors that pass through a common point. So, we 
would compute an additively weighted compact Voronoi diagram under a convex distance function instead of 
the additively weighted Voronoi diagram under a convex distance function. We will show that we can locate 
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empty (or non-empty) cells in the compact representation also. McAllister et al. [S] showed that one could 
compute a compact Voronoi diagram of the point set P having 0(|P|) complexity in 0(|P|(Iog \P\ +log |M|)) 
time. In this diagram, the closest neighbor of each query point is not a unique candidate but one of two 
candidates. To develop the concept of the compact diagram under additive weights, we need the following 
preliminaries. 




Figure 2: The spoke diagram (solid lines) depends only on the number of sites, and not the complexity of the 
distance function. The dashed lines indicate the Voronoi diagram and the solid lines indicate the compact 
diagram. 



3.1 Geometric Preliminaries 

We can give a geometric interpretation of the convex distance function between two points p, a G 
using a convex set M. Let denote the convex set M scaled by dM{p,a) and translated to p; i.e. 
M; = dM{p, a)M + p. 

Definition 5 spoke(p,a) is the line segment pa such that a lies on M^. 

Definition 6 A set X G IR^ is star shaped with respect to a if a C X and every spoke(p,a) with p G X, is 
contained in X. 

Any point k on the bisector between any two sites pi, pj G P in an additively weighted Voronoi diagram 
satisfies dM^pi, k) +uji — dM{Pj, k) where LOi and ojj are the weights of the points pi and pj respectively. 
Now, proceeding along the lines of McAllister et al. [S], we can show that the bisector between any two sites 
is a continuous curve and it separates the plane into two regions - one star shaped with respect to pi and the 
other star shaped with respect to pj . This in turn leads to the fact that a Voronoi cell of pi is star shaped 
with respect to pi. The boundary of the Voronoi cell of pi consists of portions of bisectors with other sites. 
A finite Voronoi vertex is formed by the intersection of two adjacent bisectors at a point that is equidistant 
from Pi and the other two sites defining the bisectors under additive weights. Two adjacent bisectors that 
may not intersect at a finite point is said to be Voronoi vertex at infinity. Again, as in Corollary 2.6 of 
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McAllister et al. |H], we can show that by introducing spokes from the finite and infinite Voronoi vertices 
around the boundary of the Voronoi cell of pi in Voronoi diagram of P, the cell of pi is decomposed into 
regions bounded by portions of a single p^pj -bisector. This follows from the fact that Voronoi cells are star 
shaped and the spokes lie in the corresponding Voronoi cell. Next, we define a compact Voronoi diagram 
using 0(|P|) line segments (independent of |M|), and this Voronoi diagram breaks the plane into spoke 
regions instead of Voronoi regions. In this compact Voronoi diagram, the closest neighbor of a query point 
cannot be uniquely determined, but two candidates can be determined out of which one will be the closest. 
Notice that each spoke region is a quadrilateral whose four corners are two Voronoi vertices and two sites. 
Each quadrilateral lies in the union of the Voronoi cells for the two defining sites. The compact Voronoi 
diagram is thus defined in terms of the Voronoi vertices and the spoke regions induced by the spokes. It 
follows from McAllister et al. that the number of Voronoi vertices and spokes would be 0(|P|). 

Lemma 4 The compact Voronoi diagram of \P\ sites under a convex distance function induced by a convex 
\M\-gon with additive weights has 0{\P\) Voronoi vertices and 0{\P\) spokes. 

3.2 Algorithm 

For the case with all weights zero, McAllister et al. [5] showed that one could compute a compact Voronoi 
diagram of 6(|P|) complexity in 0(|P|(log |F| + log |M|)) time such that the closest neighbor of each point 
was one of two candidates. The idea is simple if one can locate the vertices efficiently. We need to simply 
draw the spokes to each Voronoi vertex from each defining site. Figure |2] shows such a diagram and implicitly 
includes a vertex at infinity as well. Notice that the spoke diagram is composed of quadrilaterals whose four 
corners are two Voronoi vertices and two sites; each quadrilateral lies in the union of the Voronoi cells for 
the two defining sites. We now show that this diagram can also be computed with additive weights. 

Two important primitives are needed for our algorithm. 

(1) Finding the distance dM{p,q) given weighted points p and q. 

(2) Finding the Voronoi vertex for three weighted sites. 

Lemma 5 The distance dM(p,q) given weighted points p and q can be found in 0(log \M\) time. 

Proofs A binary search on the vertices of the convex polygon M finds dMip,q) and hence it requires 
0(log|M|) time. 

Next, we show the method of finding the Voronoi vertex. 

Lemma 6 A Voronoi vertex for three weighted point sites can be computed in 0(log|M|) time under a 
convex distance function du- 

Proof 6 Let P, Q and R be the sites whose corresponding weights are ojp^ojQ and ujn respectively. For 
finding the vertex of the weighted voronoi diagram for a set of sites, we consider a set of circles whose 
centers are the set of sites and radii are the weights corresponding to the sites. Assume that the sites P , Q 
and R are ordered in clockwise direction and P^, Rc and Qc be the respective circles. We want to compute 
a vertex v such that the smallest homothet of M centered at v contains the circles Pc,Qc and Rc. We draw 
the common inner tangents in clockwise direction from Pc to Qc and Qc to Rc. If there is an arc of Qc lying 
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between these two tangents, the vertex of these sites will be at infinity. Othewise there will be finite vertex. 
Now we only deal with those portions of M that touch the circles Pc and Qc when the smallest homothet of M 
centered at v contain these circles. Typically, this contact point will be an edge e of M. Compute the outer 
tangents to Pc (respectively Qc) that are parallel to e's neighboring edges. The clockwise circular arc between 
the tangential points are the probable portion of Pc that touches the smallest homothet of M containing Pc, 
Qc and Rc. Using the tentative prune and search technique of Kirkpatrick and Snoeyink J^, we can compute 
the fixed point and the voronoi vertex in 0(log \M\) time. 

Our randomized incremental construction of the compact Voronoi diagram follows the same randomized 
technique as of McAllister et al. It maintains a conflict history DAG where the nodes of the DAG 

correspond to the spoke regions. We initialize with the Voronoi diagram of three sites. To insert a new site 
p into the spoke diagram for A: < |M| sites, we locate the quadrilateral shaped spoke region containing p. 
This corresponds to two sites pi and pj any one of which may be the closest site. We measure the distances 
dAiiPijP) and dM{Pj,p) to the two existing sites that define the quadrilateral, and find out whether 
p has a non-empty Voronoi cell. If the cell is non-empty, then the Voronoi cell for new site p will carve out 
a tree from the existing Voronoi diagram; the tree topology still shows which Voronoi vertices are together 
in quadrilaterals of the spoke diagram. We can explore this tree using a number of distance and Voronoi 
vertex computations that is proportional to the number of spokes that are added or deleted. 

If we randomly order the sites and then perform incremental construction, the algorithm locates each site and 
constructs an expected 0(|P|) spokes. Thus, the total expected time is 0(|P|(log \P\ + log |M|)). Coupled 
with the initial computation of convex hull of S, and the fact that |P| = up and |5| = ns, we have the final 
result. 

Theorem 7 The set of non- dominated or skyline points of a set P of np points with respect to a set of sites 
S of ns points can be found in 0{{ns + np)logns + nplognp) expected time. 

4 Conclusions 

In this paper, we proposed an algorithm for finding the non-dominated points among a point set P with 
respect to a set of sites S in IRl^. This problem was initially proposed by Sharifzadeh and Shahabi [TT] and 
termed as spatial skyline query problem. We give some geometric insights into this problem to design an 
efficient algorithm. It would be worthwhile to extend the algorithm to higher dimensions. We intend to work 
on the dynamic version of the problem where the set of skyline points changes dynamically due to insertion 
and deletion of sites and data points. Finding non-dominated points under different domination relations 
between points will also be interesting to investigate. 

References 

[1] Stephan Borzsonyi, Donald Kossmann, and Konrad Stocker. The skyline operator. In Proceedings of 
the 17th International Conference on Data Engineering, pages 421-430, Washington, DC, USA, 2001. 
IEEE Computer Society. 

[2] K. Q. Brown. Geometric transforms for fast geometric algorithms. Ph.D. thesis. Dept. Comput. Sci., 
Carnegie-Mellon Univ., Pittsburgh, PA, 1980. Report CMU-CS-80-101. 



10 



[3] J. Cassels. An Introduction to the Geometry of Numbers. Springer- Verlag, 1959. 

[4] Jan Chomicki, Parke Godfrey, Jarek Gryz, and Dongming Liang. Skyline with presorting. In Proceedings 
of the 17th International Conference on Data Engineering, pages 717-816, Washington, DC, USA, 2003. 
IEEE Computer Society. 

[5] D. Kirkpatrick and J. Snoeyink. Tentative prune-and-search for computing fixed-points with applications 
to geometric computation. Fundam. Inform., 22:353-370, 1995. 

[6] Rolf Klein. Concrete and Abstract Voronoi Diagrams, volume 400 of Lecture Notes in Computer Science. 
Springer, 1989. 

[7] Donald Kossmann, Frank Ramsak, and SteS'en Rost. Shooting stars in the sky: An online algorithm 
for skyline queries. In Proceedings of VLDB, pages 275-286, 2002. 

[8] M. McAllister, D. Kirkpatrick, and J. Snoeyink. A compact piecewise-linear Voronoi diagram for convex 
sites in the plane. Discrete Com,put. Geom,., 15:73-105, 1996. 

[9] Dimitris Papadias, Yufei Tao, Greg Fu, and Bernhard Seeger. Progressive skyline computation in 
database systems. ACM Transaction on Database System, 30(l):41-82, 2005. 

[10] J.-R. Sack and J. Urrutia. Handbook of computational geometry. North-Holland Publishing Co., Ams- 
terdam, The Netherlands, The Netherlands, 2000. 

[11] Mehdi Sharifzadeh and Cyrus Shahabi. The spatial skyline queries. In VLDB '06: Proceedings of the 
32nd international conference on Very large data bases, pages 751-762. VLDB Endowment, 2006. 

[12] Wanbin Son, Mu-Woong Lcc, Hcc-Kap Ahn, and Scung won Hwang. Spatial skyline queries: An efficient 
geometric algorithm. CoRR, abs/0903.3072, 2009. 

[13] Kian-Lee Tan, Pin-Kwang Eng, and Beng Chin Ooi. Efficient progressive skyline computation. In VLDB 
'01: Proceedings of the 27th International Conference on Very Large Data Bases, pages 301-310, San 
Francisco, CA, USA, 2001. Morgan Kaufmann Publishers Inc. 



11 



